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Narrow-band invisibility in an otherwise opaque medium has been achieved by electromagnetically 
induced transparency (EIT) in atomic systems . The quantum EIT behaviour can be classically 
mimicked by specially engineered metamaterials via carefully controlled interference with a ’’dark 
mode” . However, the narrow transparency window limits the potential applications that require a 
tunable wide-band transparent performance. Here, we present a macroscopic quantum superconduct¬ 
ing metamaterial with manipulative self-induced broadband transparency due to a qualitatively novel 
nonlinear mechanism that is different from conventional EIT or its classical analogs. A near complete 
disappearance of resonant absorption under a range of applied rf flux is observed experimentally and 
explained theoretically. The transparency comes from the intrinsic bi-stability of the meta-atoms and 
can be tuned on/ off easily by altering rf and dc magnetic fields, temperature and history. Hysteretic 
in situ 100% tunability of transparency paves the way for auto-cloaking metamaterials, intensity 
dependent filters, and fast-tunable power limiters. 

I. INTRODUCTION 

Controllable transparency in an originally opaque 
medium has been an actively studied topic. Among those 
efforts, electromagnetically induced transparency (EIT) in 
three-level atomic systems is one of the most compelling 
ideas m- A classical analog of EIT atomic systems is 
realized in structured metamaterials that have a radiative 
mode coupled to a trapped mode 120 . The interaction be¬ 
tween two modes induces a narrow transparency window 
in which light propagates with low absorption, also cre¬ 
ating strong dispersion and a substantial slowing of light 
W- Superconducting metamaterials, especially, exhibit 
intrinsic nonlinearity thus opening the door for tunable 
electromagnetic transparency [MU- Nevertheless, most 
techniques that induce transparency require careful ma¬ 
nipulation of coupling, and only allow light to propagate 
with perfect transmission in a narrow frequency band¬ 
width inside a broad resonance feature. The narrow-band 
invisibility prevents conventionally induced transparency 
from applications requiring broad-band invisibility. In 
this work, we demonstrate a microwave metamaterial with 
self-induced broadband and tunable transparency that 
arises from an altogether new mechanism. 

Our macroscopic quantum superconducting metamate¬ 
rial is made of Radio Frequency Superconducting Quan¬ 
tum Interference Device (rf-SQUID) meta-atoms. An 
rf-SQUID is a macroscopic quantum version of the split 
ring resonator in that the gap capacitance is replaced 
with a Josephson junction. The rf-SQUID is sensitive to 
the applied rf and dc magnetic flux, and the scale of this 
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response is the flux quantum d>o = h/2e = 2.07 x 10~ 15 
Tm 2 , where h is Planck’s constant and e is the elemen¬ 
tary charge. The rf-SQUID combines two macroscopic 
quantum phenomena: magnetic flux quantization and the 
Josephson effect m, making it extremely nonlinear and 
tunable [HUMUS]. Incorporating rf-SQUIDs into meta¬ 
materials has received increasing attention [MI]. The 
large-range resonance tunability by dc magnetic flux at 
low drive amplitude [IBM], and switchable multistability 
at high drive amplitude m of rf-SQUID metamaterials 
have been studied experimentally, however, the broadband 
response of the rf-SQUID metamaterial under intermedi¬ 
ate rf flux range has not been systematically examined 
yet. 

Here we demonstrate that in the intermediate rf flux 
range the rf-SQUID metamaterial develops self-induced 
broadband and tunable transparency that arises from the 
intrinsic nonlinearity of rf-SQUIDs. Both experiment and 
simulation show that the resonance of this metamaterial 
largely disappears when illuminated with electromagnetic 
waves of certain power ranges. We can adroitly control the 
metamaterial to be transparent or opaque depending on 
stimulus scanning directions. The degree of transparency 
can be tuned by temperature, rf and dc magnetic field. We 
also discuss analytical and numerical models that reveal 
the origin of the effect and how to systematically control 
the transparency regime. The nonlinear transparency 
behaviour should also extend to the quantum regime of 
superconducting quantum metamaterials interacting with 
a small number of photons 0 [221425] . The observed tun¬ 
able transparency of the rf SQUID metamaterial offers a 
range of previously unattainable functionalities because 
it acts effectively as a three-terminal device. New applica¬ 
tions include wide-band power limiters for direct-digitizing 
rf receivers 120, gain-modulated antennas rf pulse 
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FIG. 2. Measured transmission as a function of frequency 
and applied rf flux for an 11x11 array rf-SQUID metamaterial, 
when the input rf flux (a) keeps increasing, and (b) continu¬ 
ously decreases, (c) Extracted Transparency values for sweep 
directions (a) and (b). The arrows denote the rf flux scanning 
directions. The temperature is 4.6 K. 


FIG. 1. (a) Schematic of the experimental setup for measuring 
the transmission of rf-SQUID metamaterial oriented perpen¬ 
dicular to rf magnetic field inside a rectangular waveguide. 
The dc magnetic field is applied by sending current into a 
superconducting coil outside the waveguide. The inset of (a) 
is the Resistively and Capacitively Shunted Junction (RCSJ) 
model of an rf-SQUID. Transmission of a single meta-atom 
at 4.6 K depending on frequency and rf flux is shown in (b) 
experiment and (c) simulation. Red (dark) feature denotes 
the resonance dip. Insets plot the transparency as a function 
of input power under 4.6 K (solid blue), 5.6 K (dashed red) 
and 6.5 K (dot line green) both for (b) experiment and (c) 
simulation. 


shaping for qubit manipulation, tunable intensity-limiting 
filters, and the novel concept of a metamaterial that spon¬ 
taneously reduces its scattering cross-section to near zero 
(auto-cloaking), depending upon stimulus conditions. 


II. RESULTS 

A. Broadband transparency at a fixed applied dc 

flux 


We measure the transmission and reflection of our meta¬ 
material samples positioned in a rectangular waveguide 
so that the rf magnetic field of the propagating wave 
is perpendicular to the SQUID loop (Fig. [I] (a)) (see 
Supplemental Material) [21] ■ Measured transmission of a 
single meta-atom as a function of frequency and rf flux 
4> r / at a temperature of 4.6 K under 0 dc flux is shown in 
Fig. [l] (b). Red features denote the resonance absorption 
dips of the meta-atom. At low input rf flux, the reso¬ 
nance is strong at 19 GHz [20] • In the intermediate rf flux 
range, the resonance shifts to lower frequency and sys¬ 
tematically fades away (|f>Vi = 0 dB) through the entire 
frequency range of single-mode propagation through the 
waveguide. At an upper critical rf flux, a strong resonance 
abruptly appears at the geometrical resonance frequency 
w geo / 2n = 13.52 GHz for a single rf-SQUID. We employ 
the nonlinear dynamics of an rf-SQUID to numerically cal¬ 


culate transmission (see Supplemental Material) shown in 
Fig. 0 (c) which shows the same transparency behaviour. 


We define a normalized transparency level that quan¬ 
titatively determines the degree of resonance absorption 
compared to the low rf flux absorption (see Supplemental 
Material). High transparency indicates a weak resonance 
absorption. The extracted transparency shows a clear 
onset rf flux for transparency and an upper critical rf flux 
determining the abrupt end of transparency (insets of Fig. 
[l (b) and (c)). The transparency approaches 1.0 between 
these rf flux values. The measurements are taken at 4.6 K, 
5.6 K, and 6.5 K; at lower temperature, both experiment 
and simulation show a larger range of transparency as 
well as a higher degree of transparency. 


B. Hysteretic transparency 


Collecting rf-SQUID meta-atoms into a metamaterial 
preserves the self-induced broadband transparency perfor¬ 
mance. Fig. [2] illustrates the transmission of an 11 x 11 
rf SQUID array metamaterial (see Supplemental Material 
for parameters) with interactions among the meta-atoms. 
The metamaterial is stimulated at fixed frequency while 
the rf flux amplitude is scanned under nominally 0 ap¬ 
plied dc flux at 4.6 K. The resonance is almost invisible 
as the input rf flux increases continuously through the 
transparency range (Fig. |(a)). 

However, a reverse rf flux scan renders an opaque be¬ 
haviour (Fig. [2] (b)): the resonance is strong across all 
rf flux values. Quantitatively, the transparency value 
reaches 0.9 in the forward sweep and is below 0.3 for 
the backward sweep (Fig. [ 2 ] (c)). We did numerical sim¬ 
ulations on this metamaterial and they show the same 
hysteretic transparent/opaque behaviour. Similar hystere¬ 
sis is also observed for measurements and simulations of 
a single rf-SQUID meta-atom. These observations mean 
that transparency can be turned on and off depending on 
the stimulus scan direction and metamaterial history. 





























3 


C. Transparency and bi-stability 


The origin of the nonlinear transparency is the intrinsic 
bi-stability of the rf-SQUID. The gauge invariant phase dif¬ 
ference of the macroscopic quantum wavefunction across 
the Josephson junction, 8(t), and its time dependence, de¬ 
termine essentially all properties of the rf SQUID and the 
associated metamaterial (see Supplemental Material). In 
simulation we found that S(t) is almost purely sinusoidal; 
the amplitudes of the higher harmonics are less than 
1% of the fundamental resonance amplitude for all drive 
amplitudes considered here. The amplitude of the gauge- 
invariant phase oscillation on resonance as a function of rf 
flux for a forward stimulus sweep ( Slh ) is lower than the 
amplitude for a reverse sweep (Shl) above the onset of 
bi-stability (Fig. [3] (a)). The lower gauge-invariant phase 
amplitude results in a smaller magnetic susceptibility and 
thus a reduction of resonant absorption. The relation 
between the resonance strength (degree of transparency) 
and Slh is shown in Fig. [3] (b) for the transparent state. 
The onset rf flux of transparency coincides with the abrupt 
reduction of 5LH~®rf slope and the onset of bi-stability. 

We can apply the Duffing oscillator approximation to 
analytically predict the onset of bi-stability. For interme¬ 
diate drive amplitude an rf-SQUID can be treated as a 
Duffing Oscillator (Kerr Oscillator) [55], which is a model 
widely adopted for studying Josephson parametric ampli¬ 
fiers (58j [55]. This approximation suggests that when the 
drive amplitude reaches a critical value, the amplitude 
of the 6(t) oscillation as a function of frequency is a fold- 
over resonance (Fig. [3] (c)), creating bi-stable oscillating 
states (see Supplemental Material). The amplitudes of 
the gauge-invariant phase difference oscillation for the 
transparent state (Str) and the opaque state (Sop) are 
calculated analytically for each rf flux, and compared to 
the amplitudes Slh and Shl calculated numerically (Fig. 
[3] (a)). Both the onset of bi-stability and the amplitudes 
of the two states match very well. The bi-stability of 
S(t) amplitudes explains the bi-stability of transmission 
observed in experiment and simulation (Fig. [ 3 ] (d)). The 
very good agreement between the Duffing oscillator and 
the full nonlinear numerical simulation allows us to study 
analytically how to enhance the transparency values and 
the transparency range. 

The onset rf flux value for transparency depends on 
several parameters. Higher resistance in the junction, 
higher capacitance, and higher critical current all give a 
lower onset rf flux for transparency (Supplemental Mate¬ 
rial). Operating the metamaterial at a lower temperature 
increases the resistance and the critical current, thus de¬ 
creasing the onset, explaining the modulation of onset 
by temperature observed in experiment and simulation 
(insets of Fig. [l|. The applied dc flux has a more modest 
effect on the onset rf flux. With a dc flux of a quarter flux 
quantum, the single rf-SQUID has an onset that is 13% 
smaller than the 0-flux case (see Supplemental Material). 



FIG. 3. Simulation results for a single rf-SQUID in the bi¬ 
stable regime. The pink arrows denote the signal (frequency 
or rf flux amplitude) scanning directions, (a) The amplitude 
of gauge invariant phase S(t ) oscillation on resonance as a 
function of driving rf flux amplitude. Red curves denote the 
full nonlinear numerical results. The Slh is the amplitude 
of S(t) when the frequency or rf flux scans from low value 
to high value, and Shl is the amplitude for a reverse scan. 
Blue curves are calculated analytically with the Duffing oscil¬ 
lator approximation. The Sop and Str denote the analytic 
amplitudes for the opaque state and the transparent state re¬ 
spectively. The gray arrow points out the onset of bi-stability, 
(b) The simulated transmission (blue curve, left y-axis) and 
the amplitude of S(t) (red curve, right y-axis) on resonance 
as a function of rf flux for the transparent case. The red 
dashed line shows the 5(t) amplitude interpolated between the 
low-power and high-power limits. The gray arrow shows the 
onset rf flux for transparency. It is the same as the bi-stability 
onset, (c) The fold-over resonance of the amplitude of 5(t) 
calculated in the Duffing oscillator approximation results in a 
bistable oscillation on resonance with amplitudes Str and Sop 
respectively at an rf flux of 0.0061?o. The dashed curve plots 
the unstable solutions in the fold-over resonance, (d) The 
numerically calculated transmission as a function of rf driving 
frequency at an rf flux of 0.006 < 3?o shows that the frequency 
range of bistability is the same as analytically predicted in 
Duffing oscillator approximation in (c). 


D. Modulation of transparency by dc flux 

The dc flux has a strong modulation of the transparency 
upper critical rf flux. Above the upper critical rf flux, the 
rf-SQUID experiences phase slips on each rf cycle and 
shows strong resonant absorption at the geometrical fre¬ 
quency uj geo /2TT [5Tj. We can determine the transparency 
upper critical edge when the driving frequency is fixed 
at the geometrical frequency co geo /2n, while rf flux ampli¬ 
tude scans from 0.001$ o to 0.1<I>o. The sudden decrease 
of transmission denotes the upper critical rf flux in dif¬ 
fering amounts of dc flux (Fig. El (a)). The numerical 
simulation is depicted in Fig. [4] (b). There is a tunability 
of over a factor of 3.5 in transparency upper critical flux 
by varying dc flux through the sample. Note that the 
entire dc magnetic field variation in Fig. [4] is only 10 nT. 
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FIG. 4. The transmission as a function of rf flux and dc flux 
at the geometrical frequency uj geo /2n of (a) experiment and 
(b) simulation for a single rf-SQUID at 4.6 K. The strong 
resonance absorption (red region) at the geometrical frequency 
determines the upper critical rf flux of transparency. The 
edge between the red region and the white region denotes 
the tuning of upper critical rf flux by applied dc flux. The 
experiment is taken with pulsed rf measurement with a duty 
cycle of 1% (see Supplemental Material). The pulse width is 2 
[is which is long enough to ensures that the rf-SQUID achieves 
a steady state. 

The result shows that at a fixed frequency (jj geo /2n the 
meta-atom can be transparent or opaque depending very 
sensitively on the rf flux and dc flux. 

Also, at a fixed rf flux near the upper critical edge 
the sample can be resonantly absorbing at u} geo or be 
transparent in the broadband frequency window depend¬ 
ing sensitively on the applied dc flux. Fig. [5] plots the 
experimental transmission as a function of dc flux and 
frequency when our sample is illuminated with an rf flux 
amplitude of 0.018<l>o and 0.032<J> o respectively. Strong 
absorption at the geometrical resonance appears around 
4>o/4 and 34>o/4 dc flux values, while being broad-band 
transparent near 0 and 0.5 dv Fig. [5] (c) shows that 
higher rf flux values pushes the rf-SQUID to be opaque 
under a larger range of dc flux, but the maximum value 
of transparency is higher than the lower rf flux case (Fig. 
[5] (c)). As the temperature increases, the transparency is 
weaker, as seen in Fig. [ 5 ] (d). The tuning of transparency 
by dc flux at a fixed rf flux indicates again a switchable 
on/off transparency behaviour with small variations of dc 
flux for our meta-atom. 


III. DISCUSSION 

Up to this point we have mainly discussed the results for 
a single rf-SQUID, because the transparency behaviour of 
the rf-SQUID metamaterial arises from the bi-stability of 
single meta-atoms. Disorder in the rf or dc flux can affect 
the degree of transparency in an rf-SQUID metamaterial 
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FIG. 5. The measured transmission at 4.6 K of a single rf- 
SQUID meta-atom as a function of dc flux and frequency 
at different drive amplitudes: (a) rf flux= 0.0184>o, (b) rf 
flux= 0.032<f>o- Dark red is the perfect transmission back¬ 
ground, light red and blue features denote the weak and 
strong resonance absorption, respectively. The extracted trans¬ 
parency values for rf flux amplitude of 0.018‘Fo (red dashed 
line) and 0.032<f>o (blue solid line) are shown for (c) 4.6 K 
and (d) 6.5 K. They both show a clear switchable transparent 
behaviour by dc flux. 

but the effect is quite small. In experiments on an 11 x 11 
array metamaterial, an intentionally introduced gradient 
of 0.11$ 0 applied dc flux across the array changes the peak 
transparency value from 0.94 (for the uniform applied flux 
case) to 0.91. This means the transparency is quite robust 
against noise and disorder. However, simulation shows 
that increased coupling between meta-atoms reduces the 
transparency range [30]. 

The observed tunable transparency of the rf SQUID 
metamaterial offers a range of new previously unattainable 
functionalities, enabling new applications. The root cause 
of the tunable transparency comes from the intrinsic 
bi-stability of rf-SQUIDs, which can be controlled in a 
number of ways. Some of these applications capitalize on 
the effective three-terminal device configuration in which 
transmission is modulated via a parameter, namely the 
rf or the dc flux channel. 

For example, one can design an rf transmitter where 
the input rf bias signal modulates the main channel rf 
transmission. Such a transmitter can be used in quan¬ 
tum computing to shape rf pulses for qubit control [3T], 
particularly through the scanning direction hysteresis. 

The self-induced transparency modulation feature can 
also lead to the presently unattainable gain-control ca¬ 
pability in SQUID arrays, which are the basis for highly 
sensitive quantum antennas m- This would require engi¬ 
neering a linear transparency/rf power transfer function 
in the rf SQUID metamaterial. 

Another application is in wideband directly-digitizing 
rf receivers, in which the rf SQUID metamaterial acts 
as an input power limiter to protect sensitive supercon- 
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ducting wideband digital-rf receivers from strong jammers 
[26]. With frequency selectivity, one can devise tunable 
intensity limiting filters to eliminate strong interferors 
while remaining transparent for other frequencies. 

The remarkable sharpness of transparency modulation 
by dc flux at a fixed applied rf flux amplitude (Fig. [ 5 ]) 
gives an opportunity to use single flux quantum (SFQ) 
logic to achieve fast transparency switching of the rf- 
SQUID metamaterial [20 E2]- This will enable a range 
of applications starting from SFQ-modulated digital com¬ 
munication transmitters to energy-efficient wireless data 
links between low-power cryogenic SFQ electronics and 
room-temperature semiconductor modules. Both of these 
applications are difficult to solve by means of conventional 
low-dissipation superconducting electronics. 

IV. CONCLUSION 

In summary, we show that self-induced broad-band 
transparency is observed in single rf SQUID meta-atoms 
and in rf SQUID metamaterials. The transparency arises 


from the bi-stability of individual SQUIDs. The trans¬ 
parency is hysteretic and switchable, and can be tuned 
by temperature, dc and rf magnetic flux, endowing the 
metamaterial with an ’’electromagnetic memory”. The 
transparency range and level can be enhanced through 
numerous parameters under experimental control. 
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Supplemental Material for Tunable Broadband Transparency of Macroscopic 
Quantum Superconducting Metamaterials 

DETAILS OF THE EXPERIMENT 

The sample sits inside a pulsed-tube refrigerator with a base temperature of 4.6 K. The temperature is controlled 
via an electric heater connected to a Lake Shore Cryogenic Temperature Controller (Model 340). Our sample is 
positioned in a rectangular waveguide (either X, Ku, or K-band) so that the rf magnetic field of the propagating wave 
is perpendicular to the SQUID loop and couples strongly to the meta-atoms (Fig. 1 (a)) [20 . A superconducting 
coil outside the waveguide provides dc magnetic field also perpendicular to the SQUID loop. A superconducting wire 
outside the waveguide can be used to create an intentional dc flux gradient on the sample. The sample is protected 
from environmental magnetic fields via several layers of magnetic shielding around the waveguide. The transmission 
and reflection signals from our sample are amplified by a cryogenic low noise amplifier (LNF-LNC6_20A) and a room 
temperature amplifier (HP 83020A), and are measured by a network analyzer (Agilent N5342A). The dc current 
through the superconducting coil and the superconducting wire is applied by a Keithley 220 programmable current 
source. 


PARAMETERS OF THE SAMPLES 

The single rf SQUID meta-atom and the 11x11 array metamaterial were fabricated using the Hypres 0.3 fiA/fim 2 
Nb/AlCU/Nb junction process on silicon substrates, and the meta-atom has a superconducting transition temperature 
T c = 9.2 K. Two Nb films (135 nm and 300 nm thick) connected by a via and a Josephson junction make up the 
superconducting loop with geometrical inductance L. The capacitance C has two parts: the overlap between two 
layers of Nb with SiC >2 dielectric in between, and the Josephson junction intrinsic capacitance. The rf SQUIDs are 
designed to be low-noise (T = 2 ^| r < 1 where T is the temperature and I c is the critical current in the Josephson 
junction, and Lp = pKp (^) 2 >> L [33]) and non-hysteretic (/3 r / = 2 ^ Jc < 1). 

The parameters for the single meta-atom: geometrical inductance of the rf-SQUID loop L = 280 pH, the total 
capacitance C = 0.495 pF, the geometrical resonant frequency w ffeo /27r = 13.52 GHz, the resistance in the junction 
R = 1780 Ohm (4.6 K), I c = 1.15/rA, /3 r f = 0.98, rf-SQUID inner diameter 200/jm, outer diameter 800/im. The 0 dc 
flux low-power resonant frequency is at 19 GHz. 

The parameters for meta-atoms of the 11x11 array: L = 55.99 pH, C = 2.1 pF, w seo /27r = 14.85 GHz, R = 500 
Ohm, I c = 5/rA, /3 r f = 0.86, rf-SQUID inner diameter 40/rm, outer diameter 160/mr, center-center distance 170/jm. 
The 0 dc flux low-power resonant frequency is at 20.25 GHz. 


NONLINEAR DYNAMICS AND NUMERICAL SIMULATION 

We treat an rf SQUID as a Resistively and Capacitively Shunted Josephson Junction (RCSJ-model) in parallel with 
superconducting loop inductance (see Fig. 1 (a)). The macroscopic quantum gauge-invariant phase difference across the 
junction S determines the current through the junction / = J c sinJ and the voltage over the junction V = ( H/2e)d8/dt. 




The flux quantization condition in an rf-SQUID relates 5 and the total flux (fr t ot through the loop: <5 = 2irn<S>tot/&o, 
where $o = h/ 2e is the flux quantum and n is any integer. Here we take n to be 1 because any 27 t added in <5 is 
essentially an integral constant that can be ignored in the dynamics [Mi. The total flux through the loop $t ot is the 
combination of applied flux ($>dc + sinuf) and the self-induced flux required to maintain flux quantization: 


'f’tot 




V 


,1V 


(SI) 


where the term in brackets is the total current through the junction, resistance R and capacitance C in the RCSJ 


model (see Fig. 1 (a)). Substituting d> tot with <&q5/2tt and V with (h/2e)dS/dt into equation (SI) and rearranging 
terms, we arrive at the dimensionless equation: 


d 2 5 1 dS . . . 

—— + —-— + d + p r f sind = (c pdc + <Prf sm S It) 

dr 2 Q geo dr 


(S2) 


where (f>d c = 27r &dc/ and <t> r f = 27r $> r f / $0 are the applied non-dimensional dc flux and rf flux, w geo = is the 
geometrical frequency, fl = w/w geo , r = U) geo t , Q geo = and fd r f = 2nLI c /§ 0 is the coefficient determining the 

degree of nonlinearity in an rf-SQUID. 

An array of N coupled rf-SQUIDs can be described as a system of coupled nonlinear differential equations [5D| : 


uS 


1 dS 


dr 2 Q qeo dr 


+ 6 + fdrf sin 6) = <j)dc + 4>rf sin Hr 


(S3) 


where 6 is an IV-element vector describing the gauge-invariant phases of the N rf-SQUIDs, (j>dc and <f> r f are iV-element 
vectors denoting the non-dimensional dc flux and rf flux applied to each rf-SQUID respectively, A is an N x N coupling 
matrix determining the interactions between the meta-atoms [30]. 


We solve equation (S2) (or equation (S3) for an array) numerically to determine S(t) in steady state. The time- 
dependent part of S(t) is very nearly sinusoidal, independent of the drive amplitude. With this, we calculate the 
dissipated power V 2 /R and the effective permeability fi r of a meta-atom (or a metamaterial). The transmission and 
reflection through a partially filled rectangular waveguide with the rf-SQUID medium can be calculated and compared 


to experiment (Section. IV) [20] . 

Under a fixed applied dc flux, the normalized transparency value is defined as a function of 4> r f (Fig. 1 and Fig. 2): 


Transparency{(j) r f ) = 1 — 


\S21,res((t>rf)\{dB) 

\S 2hres (0)\{dB ) 


(S4) 


where |<S , 2 i,res(^’r-/)| is the transmission on resonance at a given rf flux, and |S , 2 i,res( 0 )| is the transmission on resonance 
when the drive amplitude is low (<f> r //<f>o < 0.001 for the single rf-SQUID meta-atom). 

When our sample is driven by a fixed rf flux amplitude, the normalized transparency value is a function of the 
varying dc flux (Fig. 5): 


Transparency ((j>d c ) = 1 — 


\S2l,res(4>dc)\(dB) 

min(\S 21 ,res((l>dc)\)(dB) 


(S5) 


where \S 2 i J res{(t > dc)\ is the transmission on resonance at each dc flux, and rnin(\S 2 i,res{4 > dc)\) is the minimum transmission 
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on resonance as the dc flux varies. 


DETAILS OF DUFFING OSCILLATOR APPROXIMATION 


Given that sin d ~ 6 — |y for small argument, we get the equation of a single rf-SQUID in the Duffing oscillator 
approximation, which can be solved analytically: 


d 2 5 

dr 2 


1 d6 

Qgeo dT 


+ (1 + /3 r f)5 



<t>dc + <j>rf sin Qt 


(S6) 


In the Duffing oscillator approximation, the natural frequency of the oscillating S when <pd c = 0 is wo = y/1 + Prf^geo', 
it determines the resonant frequency of an rf-SQUID when the drive amplitude is very low. For the single meta-atom, 
wo/27t = 19 GHz. The /3 r f/6 prefactor of the nonlinear term <5 3 will modify the resonance properties when the drive 
amplitude is high. 


We separate the dc and ac parts of 6(t) for analysis. With the ansatz <5 
coupled equations for Sdc and 5 r f. 


&d,c + $ r f in equation (S6), we get two 


— ^~Q~^dc + (1 + /3rf)3dc - 


= 0 


(S7) 


d 2 5 r f 

~d^~ 


dS r f 


d T + (1 + Prf — '-^-5d c )S r f 


Qg 


- ~ ~W~ SdcS rf = 't’rf sin Hr 


Prf 


(S8) 


The non-zero (j>dc gives rise to a non-zero 5d c , which adds a S 2 ^ term in the ac gauge invariant phase oscillation. Also, 
the resonant frequency at low drive with a non-zero dc flux is modified to = uj geo y/l + (3 r f — (/3 r f/2)S 2 c (smaller 
than wo)- 

If we make the ansatz that S r f = b sin (f2r + a) [35], one finds a cubic equation for the amplitude b 2 , 

+ = (so) 

^Wgeo 

where e = w — wi is the difference between driving frequency and the natural frequency, and k is the anharmonicity 
coefficient given by: 


Prf^geo f , 5 PrfU^eJdc s 

16wi 3 uif 


(S10) 


The number of real roots of equation (S9) changes with driving frequency and driving amplitude. When the driving 
amplitude is very small, equation (S9) has one real root for b 2 throughout the whole frequency range. The peak value 


of b denotes resonance. As driving amplitude of the meta-atom increases, the resonance bends towards the lower 
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frequency side, but is still single-valued. After the driving amplitude reaches a critical value (f> r f,bi determined by 


&rf,bi 


= 2 , 


wr 


3-\/3|ac| Q^ geo oj geo 


(Sll) 


equation (S9) has three real roots for b 2 in a range of frequencies (Fig. 3 (d)). The middle-value root is unstable, so 
finally we get bi-stability in a certain range of frequencies at a high enough driving amplitude. 

With a non-zero dc flux, the oscillating S now has a dc part, and the natural frequency decreases to oj i, which 
further modifies the anharmonicity coefficient k. Since the onset of transparency predicted by Eq. |S 1 1| is related the 
resonant frequency oji and k, we can tune the onset of transparency by dc flux. For example, with a dc flux of a 
quarter flux quantum, the single-SQUID has an onset 4> r f,bi/ 27r that is 13% smaller than the 0-flux case. 


EFFECTIVE PERMEABILITY 


The magnetic susceptibility of the rf-SQUID metamaterial can be calculated as < < f > a c/‘fV/ sin fir >, where <f> ac is 
the time dependent part of the flux induced in the SQUID. Further we can calculate effective magnetic permeability 
by knowing the filling factor F of the metamaterial inside a waveguide. 


Hr — 1 + F 




sin Dr 


- 1 


(S12) 


We calculate Hr as a function of frequency as the drive amplitude through the rf-SQUID metamaterial varies (Fig. 
|S1| ). In the transparent regime the real part of Hr approaches 1.0 and the imaginary part reaches 0.0, indicating 
that the metamaterial can be manipulated by rf flux to become magnetically inert. The steps in $ r //$o = 0.3 case 
indicate the multi-stability at high drive amplitude predicted in previous work |21j . For a filling factor of only 8 %, the 
metamaterial already shows a large negative effective Hr- 

The transmission through a partially filled rectangular waveguide can be calculated as 


So 1 — 


V7 


cos kl — 


(i+ 7 ) 


sin kl 


(S13) 


where l is the effective length of the medium, k = \JHr ( 7 ) ~ ~ (a)~ the wave number in the medium, fc 0 = 
\J (y ) 2 — (f ) 2 is the wave number in the empty waveguide, c is the speed of light, a is the waveguide dimension, 
and 7 = k ko . Conversely, the effective Hr can be extracted from measured complex transmission data by employing 


equation (S13l. The Hr extracted from experimental data shows the magnetically inert behavior in transparent regime 
as well. 


TRANSPARENCY RANGE TUNING 


We evaluate the onset < fVj/ < F 0 of transparency by the Duffing oscillator approximation, and the upper critical 
<F r //<I?o by methods discussed in [21]. The two limits of transparency can be tuned by different parameters (Fig. S2|. 
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Frequency (GHz) 



FIG. SI. The (a) real part Real[/j, r ] and (b) imaginary part Imag[g r \ of the effective magnetic permeability g r of an rf-SQUID 
metamaterial in simulation. The effective g r as a function of frequency is plotted for input rf flux = 3 x 10“ 4 (Blue 

line), ‘Fr.f/'l’o = 0.132 (Red line), < £> r // l ho = 0.3 (Green line). The parameters: L = 280 pH, C = 0.49 pF, R = 1780 Ohm, 
F = 8%, rf-SQUID inner diameter: 100 (im, outer diameter: 800 pm. 



FIG. S2. The transparency rf flux range tuned by capacitance (green line, R = 1780 ohm, I c = 1.15 pA), resistance in the 
junction (black solid line, C = 0.49 pF, 7 C = 1.15 pA), and the critical current (red dashed line, R = 1780 ohm, C = 0.49 pF). 
In all cases L = 280 pH, rf-SQUID inner diameter: 100 (im, outer diameter: 800 /im. 


Larger critical current I c increases significantly the upper critical o> while the higher capacitance and resistance 

reduce the onset d> r y/d>o for transparency. The transparency range greatly broadens when either of these three 
parameters is increased. At lower temperature, the critical current and resistance increase simultaneously, resulting in 
a wider transparency o range. The enhancement of transparency range at lower temperature is consistent with 

both experimental and numerical results. 


BI-STABILITY OF FOLD-OVER RESONANCE 


In the bi-stability regime, a frequency scan and rf flux amplitude scan both result in a hysteretic transparency 
behavior. The fold-over resonance (Fig. [S3] (a)) predicted in the Duffing oscillator approximation shows that if the 
driving frequency scans from low to high, the amplitude of 5 follows the trace F-E-D-B-A, where the largest amplitude 
is at point B, denoting the transparent resonance state ( Str )■ A reversed frequency scan of trace A-B-C-E-F makes 
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FIG. S3. The analytically calculated fold-over resonance of an rf-SQUID meta-atom in the Duffing approximation, (a) The 
fold-over resonance when 4 > r //‘f > o — 0.006. The region between B-C-E-D is hysteretic. Depending on scanning direction, the 
resonance can either happen at point B with amplitude of Str or at point C with amplitude of Sop- (b) The calculated fold-over 
resonances for two rf flux values: 4 > r // < I > o = 0.003 (red line), and < h r //4>o = 0.006 (black line). Dashed lines denote the unstable 
state. Ai ( Bi ) and A 2 (B 2 ) are the lower-amplitude roots in the bi-stability regime for 4> r //4 >0 = 0.003 ( < f? T .//5 > o = 0.006) case. 
A' 2 is the higher-amplitude root in the bi-stability regime for < f? r //<f?o = 0.003 at the same frequency of A 2 . B[ and B 2 are the 
higher-amplitude roots in the bi-stability regime for <h r //$o = 0.006. 


point C the dissipative resonance (Sqp)- 


The hysteresis in rf flux is similar. Fig. S3 (b) plots the bi-stable fold-over resonance of an rf-SQUID meta-atom 
at drive amplitude & r f /&0 = 0.003 (red curve) and d> r //d>o = 0.006 (black curve). As the rf flux amplitude keeps 
increasing at a fixed frequency, in this case 17.6 GHz, the amplitude of 5 changes from Ai to Bi, and stays in the 
transparent state. However, if the scan starts from very high rf flux value and gradually decreases, the higher amplitude 
B( which denotes the dissipative mode will be excited. Further reduced rf flux $,.//$0 = 0.003 brings the rf-SQUID 
back to the A\ state. The highest amplitude at a fixed frequency (17.6 GHz) is at B\ when sweeping rf flux from low 
to high values, while it is at B( for a reverse scan. The same applies for the rf flux amplitude scan at 18.5 GHz: a 
forward sweep in rf flux amplitude keeps the amplitude of S in the lower values (A 2 to B 2 ), while the backward sweep 
excites S to oscillate at B' 2 and A' 2 , resulting in a dissipative state. 


PULSED-RF MEASUREMENT 

In Fig. 1 in the main text, the experiment and simulation match well except that the upper critical rf flux differs 
by about 0.12d>o. We believe the difference is due to the local heat generated on our sample from high input power 
that turns the transparency state to the dissipative phase slip state prematurely. To reduce the local heating, we 
conducted a series of pulsed-rf measurements, where the input signal is periodically pulsed with the pulse width long 
enough so that S(t) achieves a steady state. Pulsed rf measurements give larger transparency and better resemblance 
to simulation. For similar reasons, the experimental data in Fig. 4 (a) in the main text was taken with the pulsed rf 
measurement with a duty cycle of 1%. 













